Breakdown of adiabaticity when loading ultra-cold atoms in optical lattices 
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Realistic simulations of current ultra-cold atoms experiments in optical lattices show that the 
ramping up of the optical lattice is significantly nonadiabatic, implying that experimentally prepared 
Mott insulators are not really in the ground state of the atomic system. The nonadiabaticity is 
even larger in the presence of a secondary quasi-periodic lattice simulating "disorder" . Alternative 
ramping schemes are suggested that improve the adiabaticity when the disorder is not too large. 
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Ultra-cold atoms in optical lattices form a wonder- 
ful toolbox ^ for creating novel matter phases realiz- 
ing a "condensed matter theorist dream" . This comes 
from the flexibility of parameters in cold atom physics 
where both the depth of the optical lattice (changing ei- 
ther the laser intensity or its detuning) as well as atom- 
atom interactions (via external magnetic field and Fes- 
hbach resonances) can be modified with unprecedented 
precision. By increasing the depth of the optical lattice 
where a Bose-Einstein condensate is loaded, one induces 
a transition between the superfluid (SF) phase and the 
Mott insulator (MI) phase, as suggested in (basing 
on the Bose- Hubbard (BH)_tight binding model [31]) and 
attempted few years later 4] . That triggered several ex- 
perimental and theoretical studies involving bosons @,[3]! 
fermions and their mixtures (see [7| for details). 

Disorder can be introduced in cold atom systems in a 
controllable way using optical potentials created by laser 
speckles or bichromatic lattices Q. This exciting possi- 
bility has attracted a lot of research For weakly (or 
non) interacting bosons, Anderson localization has been 
recently observed For strongly interacting bosons 

in a disordered lattice, studies of the BH model (1, |Tl| 
revealed the existence of a novel insulating but compress- 
ible phase called Bose glass (BG) phase. 

As far as we know, a clear demonstration of the exis- 
tence of the BG phase in traditional condensed matter 
experiments is yet to be made. In a recent experiment 
[12'] using ultra-cold atoms, a smearing of the absorption 
peaks was observed in the presence of "quasi-disorder" 
and interpreted as being consistent with the existence of 
a BG. While the ID BH model in bichromatic lattice 
has been discussed theoretically in detail [13], no direct 
comparison of the theory and experiment has been made. 
Such an analysis is therefore needed. 

The system can be described by the BH Hamiltonian: 



H 



J y: bib,, 

U-j') 



u 



where bj (fej) is the destruction (creation) operator of one 

particle in the j-th site, hj = b^jbj is the number oper- 
ator^ and (j, j') indicates the sum on nearest neighbors 
0, Q . The interaction energy U and the hopping energy 
J can be evaluated via appropriate integrals of Wannier 
functions obtained for lattices of a given depth. 

The sites energies, e^, are given by the sum of the har- 
monic trap potential and - if present - an additional 
"quasi-disorder" due to a second weak optical lattice cre- 
ated by a standing wave with wavelength A2 : 
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where m is the particle mass, a the lattice spacing, lu the 
trap frequency, jo the position of the trap center, Ef(2 
the recoil energy associated with wavelength A2 I4] , and 
(f) the phase of the A2 standing wave. S2 measures the 
amplitude of the disorder in recoil energy units 14 ] . 

In the experiment [l^ . an array of ID systems is ob- 
tained by loading a ^^Rb condensate in a 2D optical lat- 
tice ("transverse lattice"). The height of the 2D lattice, 
s±_ ~ 40, results in a negligible tunneling rate between 
the one-dimensional atomic tubes. During the ramping 
on of the "transverse lattice" , an optical potential in the 
third direction is turned on as well, along the axes of 
the atomic tubes. Thus a collection of independent ID 
atomic samples is loaded in the superposition of a "main" 
optical lattice (at Ai = 830 nm with a height si = 16) 
and a "disordered" optical lattice (at A2 — 1076 nm with 
a height S2). All lattices are slowly ramped up, increasing 
the intensity of the laser beams following an exponential 
ramp over 100 ms time. After the loading, the excitation 
spectrum of the obtained state is measured by modulat- 
ing the intensity of the "main" lattice beam for 30 ms 
with an amplitude of 30%. Also, the coherent fraction is 
extracted from time of flight density profiles 0, H, [12] . 

Our numerical simulation of the experiment is per- 
formed using a TEBD algorithm [3] (also known as t- 
DMRG algorithm It has already been applied to 
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the BH model 17|, [l^ ; we use it here for parameters of 
a given experiment aiming at reproducing the processes 
taking place during ramping up of the lattices potential. 
Since the theory is based on the BH model, it cannot sim- 
ulate the very initial stage of the experimental procedure 
that starts from a harmonically trapped Bose-Einstein 
condensate. In the first stage, the atomic gas is in a SF 
state; ramping up the lattice must be done sufficiently 
slowly to ensure that all atoms remain in the lowest band 
of the periodic potential. The optical potential must be 
increased from to one recoil energy in a time longer 
than the recoil time (300 /xs), which is largely fulfilled in 
the experiment. The atomic wave function then follows 
adiabatically the shrinking of the Wannier state in each 
site as long as tunneling between neigboring sites remains 
dominant. Following this reasoning, our simulations are 
performed starting from si = 4 or 5, where the system 
is already in the tight binding regime, in the SF phase. 
We have checked that our results are practically inde- 
pendent of this choice (proving that the initial stage is 
adiabatic). Wc choose 151 particles occupying roughly 80 
lattice sites, in order to match the experimental parame- 
ters of [l^ . Similar conclusions are obtained for different 
values of the parameters. 

In the experiment IJI, an absorption peak at fre- 
quency 1.9 kHz was observed, while the interaction en- 
ergy for si — 16 is U — 2.3 kHz. We thus believe that 
the si parameter was overestimated. To match the po- 
sition of the peak, we modify all laser intensities by a 
7/8 factor yielding si = 14 and s± = 35. Attempts to 
simulate the absorption curve for the 30 % modulation 
of the lattice lasting for 30 ms failed. Such a strong 
modulation apparently creates too much excitation and 
entanglement across the lattice for the TEBD algorithm 
to be reliable. The absorption spectra obtained for 10-30 
times smaller modulation resemble those obtained in [l3| 
with a prominent peak around U and a much smaller one 
around 2U (in experiment this relation is reversed). In 
the presence of disorder, the absorption peaks broaden — 
both in the real experiments and in our numerical simu- 
lations — but it is hard to extract any information from 
this fact: firstly because of nonlinearities induced by the 
experimental strong modulation, and secondly because 
absorption spectra turn out to be almost insensitive to 
the details of the dynamics such as a breakdown of adi- 
abaticity. 

The numerical simulations make it possible to ana- 
lyze the state of the system by different means. The 
first step in this direction is to test an assumption done 
in ultra-cold atom " quantum phase transition" exper- 
iments d, H, [13], namely the adiabaticity of the opti- 
cal lattice switching on. This assumption is especially 
critical in the laser intensity region corresponding to the 
transition from the SF to the MI phase. For the infinite 
system, due to the critical slowing down [l9| , such a tran- 
sition cannot be fully adiabatic. For large 3D systems f^. 



an approximate mean-field simulation [20| showed that 
indeed adiabaticity is broken. For a collection of ID 
tubes 0, [l3] , the small system size may assure the adia- 
baticity of the transition/crossover between phases. We 
can directly test the adiabaticity by comparing the state 
obtained dynamically by our procedure with the exact 
ground state of the system. The latter can be found by 
imaginary time propagation (using essentially the same 
algorithm). Next, the overlap between the dynamically 
created state and the ground state can be calculated. 

This procedure yields for a single lattice (s2 = 0) an 
overlap equal to 0.095; for weak disorder the correspond- 
ing overlap drops to 0.02-0.05, depending on the value 
of S2- For S2 greater than unity, the overlap drops dra- 
matically to 10""^ or smaller values (note that S2 ~ 1 
corresponds to disappearance of the MI phase accord- 
ing to BH model predictions). The overlap values ob- 
tained indicate that the switching on of the lattices is 
not adiabatic. What is the character of the state pre- 
pared at the end of the ramping up? It is some kind 
of "wavepacket" , a coherent superposition of the ground 
state and some excited states. A crude way of measuring 
the wavepacket character is to compute the "excess en- 
ergy", i.e. AE — ^^(wavepacket)— £^(ground state). Re- 
gardless of S2 it is about 0.15J7 meaning that the system 
is not highly excited. 

More precise information can be obtained by con- 
sidering the temporal autocorrelation function of the 
wavepacket. Indeed, the wavepacket can be expanded 
onto the eigenstates \(pi) (with energy Ei) of the Hamil- 
tonian at the final values of si and S2- \'^) = '^ilfi)- If 
we now let j^") evolve freely with the final Hamiltonian, 
we obtain the autocorrelation function: 



C(() = (*(0)|»(l» = i;i<.,P<!xp 



Thus Ei 's and overlaps can be obtained by Fourier trans- 
forming C(t). For finite time evolution, a state of the art 
harmonic inversion technique [2l| provides us with very 
accurate results. The exemplary results are shown in 
Fig. [TJ In the absence of disorder, about 8 states con- 
tribute significantly to the wavepacket with surprisingly 
one excited state having 30% overlap. For strong dis- 
order, the dynamically created wavepacket spreads over 
several tens of eigenstates. There are large fluctuations 
of the overlaps among neigboring states, meaning that 
the excitation brought by non adiabatic switching dif- 
fers from a thermal excitation. Note also that real-time 
TEBD propagation brings here a spectral information 
not accessible to standard "quasi-exact" numerical meth- 
ods such as DMRG or Quantum Monte-Carlo. 

Let us underline that, while we have analyzed as closely 
as possible the ramping procedure in the recent "disor- 
der" experiment |12l |. our findings should be applicable 
to other similar experiments in which both lattices and 
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FIG. 1: (color online) Spectral analysis of the wavepacket ob- 
tained dynamically for si = 14 after switching on the lattice. 
The solid line is the Fourier transform of the autocorrela- 
tion function, eq. ((3)l, showing peaks at the energy levels of 
the system, with an intensity equal to the squared overlap 
with the wavepacket (an exponential ramp is used). About 
eight (upper panel, no disorder) or several tens (lower panel, 
S2 = 0.4375) of states are significantly excited proving that 
the preparation is not adiabatic. Crosses indicate the values of 
overlap obtained via a harmonic inversion technique. Filled 
red circles are the corresponding overlaps for the improved 
ramp discussed in the text. 



the harmonic trap are rainped on time scale of a 100 ms 
either in ID Q or in 3D [illl. 

An obvious way to improve adiabaticity is to make 
the ramping much slower. That leads, however, to other 
problems (e.g. decoherence due to spontaneous emission) 
- already a 100 ms ramping time seems quite long. An- 
other possibility is to optimize the ramping profile. The 
main source of nonadiabaticity is likely to come from 
the vicinity of the SF-MI phase transition. While, for a 
homogeneous BH model, si ~ 5 corresponds to the tran- 
sition, the presence of the trap shifts the transition point 
higher in si. 



see 



2|). 



Inspection of ground states for 
various si shows that the crossover point is si !=s 8, in 
agreement with experimental observations [l^ . There- 
fore, a ramp which slows down around si = 8, for ex- 
ample si{t) = 8sin^(27rt/T), and again speeds up to- 
wards the final value (for example using an exponential 
increase) is likely to improve adiabaticity. This is fully 
confirmed by our numerical simulations shown in Fig. [2l 
The "best" pulse shape slowing down around si = 8 per- 
forms quite well, allowing to reach 53% overlap with the 
ground state, a considerable improvement over the expo- 
nential ramp. The excess energy for this pulse is more 
than twice smaller than for the exponential pulse. Other 
shapes perform almost as well provided the region around 
si = 8 is passed slowly. In contrast, the "sigmoidal" 
ramp of the form si{t) = 14/[1 + exp(— a(t — tmid)/T)] 
where T = 100ms is the total ramp time, imid = T/2, 
a = 20 as in [2^, with a fast change in the vicinity of 
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FIG. 2: (Color online) Top panel: Different ramps used for 
turning on the lattices over 100 ms. The solid curve corre- 
spond to the experimental exponential ramp with time con- 
stant T = 30 ms. The dashed blue curve is the "best" pulse 
with a sin^ ramp up to to si = 8 in 50 ms, followed by an 
exponential increase with time scale r = 10 ms. The dotted 
red curve is the sigmoidal ramp used recently [23, . Bot- 
tom panel: Overlaps of the wavepacket obtained at the end 
of the ramp (starting from the ground state at si =4), with 
ground states at different si values. The evolution remains 
reasonably adiabatic up to si = 7. The "best" pulse performs 
significantly better than the exponential ramp yielding more 
than 50% final overlap. The sigmoidal ramp is the worst one. 



the phase transition, leads to a total loss of adiabaticity. 
For experiments exploring the insulator phase 23], this 
is certainly not a good choice. 

While the current experiments exploring the MI regime 
appear not to be strictly adiabatic, this affects slightly 
properties such as occupation of different lattice sites 
or momenta distributions. The dynamically created 
wavepackets thus share with the corresponding ground 
states some global properties. In this approximate sense, 
the creation of MI state in ID experiments without dis- 
order [5] is confirmed by our simulations. 

The situation is quite different in the presence of the 
secondary lattice creating quasi-disorder. Already for a 
relatively small disorder, S2 = 0.4375, a few excited states 
are significantly populated, see Fig. [TJ The situation be- 
comes drastically worse for S2 > 1 — the region where MI 
should disappear completely. There, the overlap between 
the dynamically created wavepacket and the ground state 
becomes vanishingly small, 10~^ or less for the exponen- 
tial ramp at S2 = 2.1875. However, the experimental 
results reported in [l3 | support the claim that a behav- 
ior characteristic of BG was observed. 

The calculated momenta distributions averaged over 
different phase shifts between lattices, show a broad peak 
similar to the one for the ground state at si = 14 and 
S2 = 2.1875. The difference between the ground state 
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While competing this work we became aware of a re- 
cent work [26[ which shows that adiabaticity is much 
harder to obtain in the presence of two lattices also for 
shallow lattices (Gross-Pitaevskii regime). 
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2004-517186). 



FIG. 3: (Color online) Average occupation (upper panel) and 
its variance (lower panel) of individual sites, for si — 14 and 
S2 = 2.1875, i.e. for parameters where a Bose glass phase 
is expected. Black circles (and solid line) correspond to the 
ground state of the system, with essentially only integer occu- 
pations and small variances, together with fluctuations of the 
occupations along consecutive sites, as expected for a Bose 
glass. Red crosses (and dashed line) correspond to the dy- 
namically created wavepacket. Non- integer occupations and 
largely increased variances show significant excitation of the 
system and that its properties are different from a Bose glass. 



and the wavepacket becomes evident looking at the occu- 



pations of sites and their variances Aui — \J — (ni)^ 
as shown in Fig. [3l The ground state is characterized by 
integer occupations with relatively low number variance, 
while the dynamically excited wavepacket shows an av- 
erage variance twice larger and non-integer occupation 
of some sites. The optimized pulse, working well for no 
disorder, leads to similarly disappointing results. 

In summary, we have shown that the experimental 
preparation of a MI in recent experiment [l^ is partially 
non adiabatic. Similar conclusion should apply to earlier 
experiments 0, [1] . Still the global properties of the cre- 
ated wavepackets resemble those for the ground states (in 
particular regions of integer occupations and low number 
fluctuations appear as in the MI phase). A careful re- 
shaping of ramps used to switch on the lattices may lead 
to an adiabatic creation of MI. 

The presence of disorder amplifies nonadiabaticity. 
Our simulations indicate that the properties of the 
ground state and of dynamically created wavepackets dif- 
fer substantially indicating that a BG was not actually 
achieved in the experiment ^12] . The observables mea- 
sured in this experiment are not sufficiently sensitive to 
the state of the system to characterize unambiguously 



a BG. To that end, other tools suggested recently [25 1 
might prove useful. A simple optimization of the pulse 
shape, successful in the absence of disorder, does not 
help. We are currently working on other schemes aiming 
at dynamical creation of a BG. 
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